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Abstract 

This paper presents an application of time-frequency 
methods to characterize the dispersion of acoustic 
waves travelling in a one-dimensional periodic or dis- 
ordered lattice made up of Hclmholtz resonators con- 
nected to a cylindrical tube. These methods allow 
(1) to evaluate the velocity of the wave energy when 
the input signal is an acoustic pulse ; (2) to display 
the evolution of the spectral content of the transient 
signal ; (3) to show the role of the localized nonlinear- 
ities on the propagation .i.e the emergence of higher 
harmonics. The main result of this paper is that the 
time-frequency methods point out how the nonlinear- 
ities break the localization of the waves and/or the 
filter effects of the lattice. 

1 Introduction 

When considering propagation of acoustical pulses in 
dispersive media, the crucial question is that of the 
spreading of the signals. Indeed, the investigation 
of the dispersion relation gives a lot of information 
about the processes which play a part in the wave 
propagation. Efficient tools for characterizing the 
temporal localization of the spectral components are 
the Time-Frequency Representations (TFR) of sig- 
nals. Joint time- frequency representations combine 
time and frequency-domain analysis by displaying a 



signal as a function defined over the time frequency 
plane. These representations make the behavior of 
transient signal visible which is, if not impossible, at 
least, very difficult by using harmonic analysis. For 
waves travelling in strongly dispersive media (e.g. a 
periodic lattice) these approaches can be very useful 
since they allow the determination of the dispersion 
relations by measurement of the time of arrival of the 
different components (or frequencies) of the input sig- 
nal. When the propagating medium contains nonlin- 
earities, the TFR lead also to a detailed visualization 
of the nonlinear effects and show their consequences 
for the propagation of signals. 

In the following, the dispersion characteristics of 
sound waves propagating in a waveguide with an ar- 
ray of Hclmholtz resonators connected axially are ex- 
amined. From a general model of propagation tak- 
ing account of weakly nonlinear effects localized in 
the resonators, we present Time-Frequency pictures 
given by acoustic pulses and frequency modulated 
signals (chirps). 



2 Propagation in a one- 
dimensional discrete medium 

2.1 Lattice description 

A one-dimensional lattice made up of an infinitely 
long cylindrical waveguide (call also pipe) connected 
to an array of Hclmholtz resonators numbered by n 
(simply call resonator hereafter) is considered. The 
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resonators are connected to the pipe through a pin- 
point connection, the radius of the throat's cross sec- 
tional area s n of the n th resonator being assumed to 
be small compared to the wave length of the acoustic 
wave (y/s^/X <1). Each connection is located along 
the axis of the waveguide by its coordinate z n with 
axial spacing d n for two consecutive points as shown 
in Fig. (|IJ. 
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Figure 1: Experimental set up. 



2.2 General case 

In a section of the pipe between two consecutive con- 
nection points, the acoustic wave described by the 
pressure p(z, t) and the acoustic velocity v(z, t) is the 
solution of the wave equation : 



d 2 p(z,t) 1 d 2 p{z,t) 



dz 2 



dt 2 



= 



(1) 



where c is the sound speed in free space. 
At each connection point (Fig. fl}, the boundary 
conditions require the conservation of acoustic flow 
and continuity of acoustic pressure : 
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where v t (z n , t) is the acoustic velocity in the throat of 
the (n + l)^ n resonator. By using the Euler equation, 



Eq. © and © become 
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where p is the air density. They lead to the inhomo- 
geneous wave equation [TJ ; 
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where the right hand side term has a comb-like struc- 
ture, the teeth of which having a length related to the 
corresponding impedance jump. This forcing term 
acts as an array of secondary point-sources (scatter- 
ers) which work when they are illuminated by the 
wave travelling in the pipe. 

2.3 Propagation of a monochromatic 
wave 

2.3.1 Propagation equation 

For a monochromatic acoustic wave with a frequency 
below the cut-off frequency of the waveguide, the 
acoustic pressure p(z, t) and velocity v(z, t) along the 
waveguide are : 
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p(z,t) = p(z)e 3ul and v (z, t) — v(z)e 

where u> = kc is the angular frequency. The ampli- 
tudes p(z) and v{z) are related by an impedance rela- 
tion. At each connection between the waveguide and 
a resonator, the wave impedance (which is modified 
by a change of cross sectional shape), the acoustic 
velocity and therefore the derivative of the pressure 
dp/dz, are discontinuous functions. Using the Eq. 
|6]). the pressure p(z) is given by the solution of the 
following equation : 



d 2 p(z) 
dz 2 



k 2 p(z) = ^^{z - z n )a n p(z), (7) 



where a n = —jujps n /{SZ n ). In this relation Z n is the 
impedance of the n th resonator connected at point 
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z n seen from the guide and a n is the jump of the which links the vector pairs (A n B n ) and 

(A n -\-iB n +i) . The matrix T n+ \ has the following 



pressure derivative -fi\ z — Zn . 
2.3.2 Matrix method 

In the n + cell (z n < z < z n+ i) the pressure and 
the acoustic velocity are respectively denoted by p n 
and v n . The solution of the Eq. ([7]) can be found by 
the matrix method [IJ[2l[3] and the solution p n (z) is 
given as a linear combination of waves travelling in 
opposite direction : 

Pn (z) = A n e^- Z ^ + B n e^ k{z - z -\ (8) 

Here the coefficients A n and B n are respectively the 
amplitudes of the forward and backward waves. 
This function is connected to p n +i in the next region 
having passed through the impedance discontinuity 
at the connection point z n . The connectivity condi- 
tions at the scattering region are the continuity of the 
pressure 

Pn(z n ) =Pn+l(Zn), 

and the continuity of the mass flux 
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These junction conditions ensure a physically sound 
solution by requiring continuity all along the waveg- 
uide. Applying the matrix method, these equations 
may be reduced to a transfer matrix T> n relating the 
amplitudes of waves across the junction z n . From 
z n + e to z n+ i — e, the propagation is modeled by a 
phase matrix M. n +i 



M n+1 



e jk(z n+ i-z n ) 






e -jk(z n +l- 



The transfer matrix of propagation over the region 
z n < z < z n+ i is T n +i — M n +iD n . In this way the 
pressure in the duct is given by the following recursive 
relation : 



V n+ i = T n+1 V n where V n 
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form [T] 



Ijk, 



a jkd„ 



(10) 



The propagation through the lattice from z n to z n + m 
is then described by the relation 



W %i+iV n 



(11) 
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V n may be interpreted as the vector of the initial 
conditions (or boundary conditions) and V n + m is the 
vector of the wave amplitudes m cells further along. 
The effect of disorder is always to break some symme- 
try [1]. In the case of our lattice disorder breaks the 
periodicity. So it is easier to characterize disordered 
systems in terms of their deviations from an ideal of 
order than it is to define a perfectly disordered sys- 
tem on which some partial degree of order is to be 
imposed, and, to think about a disordered system it 
is necessary to keep in mind the ideal system from 
which it derives. So, before tackling the problem of 
wave propagation in disordered media, some specific 
results about periodic lattices are recalled. 

2.3.3 Periodic lattice 

For an ordered (i.e. periodic) lattice, the transfer ma- 
trix T n = T is the same for every cell, and the infinite 
medium problem is analogous to the Kronig-Penney 
model well known in solid state physics to investigate 
the motion of electrons in a periodic potential [5] . So 
the propagation in the lattice can be seen in terms 
of plane waves subject to multiple reflections at each 
derivation, resulting in standing waves. It is also pos- 
sible to describe the propagation in terms of a collec- 
tive excitation that propagates in the periodic lattice 
without scattering but with a modified dispersion re- 
lation. The result is that spatial periodicity gives rise 
to dispersion even in the model of plane waves and 
in this special case the spectrum shows frequency in- 
tervals (gaps) where no energy propagates. 
Studying the spectral properties of such a system is 
then to seek stationary states of excitation that sat- 
isfy prescribed conditions at one end of the lattice. 
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The general theory of eigenvalues assures that the 
spectrum is not significantly affected by the values 
chosen as boundary conditions, provided the number 
of cells is large. The matrix T has two eigenval- 
ues a ± with two eigenvectors W ± respectively. Any 
acoustic wave in the duct can be represented as a 
linear combination of these eigenvectors 

p n = u + W + + u~W~. 

The effect of operating on p n by the transfer ma- 
trix T depends on the nature i.e. real or complex of 
the eigenvalues cr^. When they are real, T simply 
" pushes the wave" in the direction of the eigenvector 
corresponding to the greatest eigenvalue. If a + > or 
then 

p n+m =T m Pn ^(a + ) m u+W+ 

as to becomes large. W + is a fixed point or an in- 
variant point of the transformation p n +\ — f(p n )- 
By iteration of this transformation, W + appears as 
an attractor (or a sink) in the subspace spanned the 
eigenvectors : almost all the waves are attracted by 
W + . When the eigenvalues are complex, the itera- 
tion produces a periodic sequence : are indiffer- 
ent points of the plane. 

Because the transfer matrix is unitary, the eigenval- 
ues of matrix T obey to the equation 

a 2 + aTr{T) + 1 = 0, (12) 

where 

TriT) = 2 cos (kd) + - sin (kd) = 2 cos(qd) (13) 
k 

represents the dispersion relation of the Bloch waves 
5 , q being the Bloch wave number. This derived 
dispersion relation exhibits the peculiar characteris- 
tic of filters marked by forbidden frequencies or gaps 
or stopbands and allowed frequencies or passbands 
in the frequency domain which results from the reso- 
nances and the periodic arrangements of the medium. 
Waves that obey the relation | cos(q<i)| < 1 are within 
a passband and travel freely in the duct and waves 
such that | cos(qd)\ > 1 are in a forbidden band and 
are quickly damped spatially. They become evanes- 
cent so that they cannot propagate. 



To deduce the particularity of the band structure, Eq. 
(|13p is written in the following form 

cosikd + 9) 
cos(qd) = — — 

cos(y) 

where 9 is defined by tan(0) = The frequencies 
which bound stopbands and passbands are given by 
the equation : 

cos(kd + 9) = ±cos(0) 

leading to two kinds of solutions kd = nix and kd = 
mr - 29 with neN. 

Two kinds of stopbands appear in the band structure 
: one kind called Bragg stopband due to the period- 
icity of the lattice, the other due to the resonances of 
the scatters and called resonance stopbands. A plot 
of the dispersion relation (fT3|l corresponding to the 
ordered lattice of Hclmholtz resonators is showed in 
Fig. [2] which illustrates the different stopbands. 




Frequency (Hz) 



Figure 2: (a) Dispersion relation (from eq. [T3|) of a 
ordered lattice, (b) Transmission coefficient of the 
lattice : — simulation (from eq. I15|) . experiment. 
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2.3.4 The case of a disordered medium 

A disordered lattice is mathematically characterized 
by a random sequence of non-identical transfer ma- 
trices T n whose overall product describes the propa- 
gation along the guide. The disordered linear system 
gives rise to the classical Anderson localization [6]. 
In this case, all the matrices are different each corre- 
sponding to a different cell, and the wave propagation 
is described by the general Eq. (fTTj) . The estimation 
of the Lyapunov exponent, from the Ossedelec theo- 
rem [7], gives the behavior of the solution p(z) in the 
limit n — > +oo : 

n 

i= Urn - In (Tr ([[%)). (14) 

i=l 

An alternative way to analyze the behavior of the 
waves is the direct evaluation of r„ >TO defined, from 
Eq. (HU, by 



corresponding to the transmission coefficient of a lat- 
tice made up of m cells. As \m\ — > oo, T n m is related 
to £ by T n m w exp (— £|m|d) where d is now the mean 
value of the set of values \zi — i.e. the mean 

length of cells [8]. When the medium is disordered, 
the gaps widen : eventually waves cannot propagate 
in such a medium. This result is well known since 
Anderson's original paper [6] : in a one-dimensional 
disordered system, almost all the elementary excita- 
tions are localized. The waves cannot propagate in 
a disordered medium because of the many scatterers 
they encounter. The possibility that a wave can be 
localized in a random medium is mysterious because 
localization involves a change in the wave character. 

2.4 Propagation in a nonlinear 
medium 

For high sound level (~ 120 dB), the relation between 
the acoustic pressure and the velocity in the neck of 
the Helmholtz resonators is no longer linear. Nonlin- 
earities result mainly from the complicated motion 
of the air near the tubes between the wave guide and 
the neck (throat) of the resonators and at the aper- 
ture of each resonator [9]. 



Consequently, Eq. © is a nonlinear wave equation 
where Vt(z, t) is not a linear function of p(z, t). Nev- 
ertheless we assume that the propagation (between 
two resonators) remains linear. So the pressure in the 
main pipe is calculated in the same way as for a disor- 
dered medium but now nonlinear operators describe 
the scattering of waves at each connection point of 
the lattice. A simple model of nonlinear Helmholtz 
resonator is developed in the following by using a 
Taylor's development of the restoring force due to 
the change of pressure in the Helmholtz cavity. 

3 A simple nonlinear 

Helmholtz resonator model 

It is well known that a simple model of the Helmholtz 
resonator requires the following assumptions : (1) the 
pressure inside the cavity is spatially uniform, (2) 
the fluid in the neck moves like a solid piston. In 
this case, the air enclosed in the resonator acts as 
a spring for the lumped mass of air moving within 
the neck. A general description of the nonlinear be- 
havior of the Helmholtz resonator may be derived by 
taking into account the quadratic term in the restor- 
ing force of the spring. For thermodynamical pro- 
cesses that occur in the air within the cavity, the 
adiabatic changes of the pressure p and the volume 
V are related by po/p = (V/Vq) 1 where the subscript 
stands for the unperturbed reference values and 7 
is the specific heat ratio. The relative change of the 
pressure p n in the cavity of the n th resonator due to 
a small displacement x n of the air in the neck induces 
a restoring force F n that has the following form [TO] 

2 2 

pC S^i r 2 /3m 

f n — Pn s n — 77 [%n Q n X n T" 0[X n )\. 

In this equation a n = (7 + l)s n /(2Vo) and c is the 
sound velocity given by c = y/ jp / p. The spring 
force is no longer linear and its stiffness is now de- 
scribed by two constants. For a monochromatic wave 
the displacement x n of the air in the neck is related 
with the acoustic velocity v n — v{z n ) by the relation 
v n = jujx n and the Euler relation applied to the air 
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mass m — pl' c s n (where l' c is effective neck length) 
submitted to the harmonic force p n /(pl' c )e^ UJt gives : 

jujv n +wg a n { — Y + o«) = — 

where = SnC 2 /(Vb^) is the resonance frequency 
of the Hclmholtz resonator. The nonlinear relation 
between the acoustic pressure and the velocity just 
outside of the opening of the n th resonator is 



Pn 



z 
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Here Z^ L is the "nonlinear impedance" of the n th 
resonator : 
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and Z% is the linear impedance of the n th resonator 

2 

(Z^ = ju>pl' c (l — ^|)). So, the nonlinear effects 
lead merely to a additive correction to the linear 
impedance of the Helmholtz resonator which is non- 
vanishing only around the Hclmholtz resonance luq. 
Others nonlinear phenomenon can occur due to the 
presence of turbulences (vortex) which are generated 
around the edges of the neck [9]. According to the 
shape (sharp or rounded) of these edges, the physi- 
cal phenomena lead to nonlinear terms more or less 
important 

4 Signal analysis 

Time-frequency representations have several advan- 
tages over conventional harmonic Fourier analysis. 
Among these one can cite their ability to analyze 
multicomponcnt signals i.e. broadband signals and 
also signals containing several modes which interfere 
or combine each other. But the most important ad- 
vantage is their capacity to track the time evolution 
of transient signals. This is exemplified by the visual- 
ization of reflected signals due to scattering waves and 
of the spreading of acoustic pulses in the joint time 
frequency plane. In this paper the quadratic TFR 
(QTFR) are used to determine empirically the disper- 
sion of the waves since this provides a direct access to 



the group velocity and allows to extract pertinent in- 
formation associated with the waves travelling in the 
lattice. Unfortunately, QTFR also have some draw- 
backs. As it is bilinear in the signal, the QTFR of 
the sum of two signals is not the sum of the QTFR of 
each of them but has an additional cross term which 
gives rise to artefacts and confuses the picture in the 
time-frequency plane. Several solutions are available 
for reducing these terms, for example, by using an 
appropriate smoothing window function and/or the 
analytic signal. 

For the complex- valued time signal x(t), the general 
form of the QTFR is [12] 



-j'Bt—jruj—j'du 



i 

27T. 

x(t+ ij)x*{t- ^jdrduM (16) 



where r) is the kernel of the QTFR, * stands 
for complex conjugate and t and ui are time and fre- 
quency respectively. The most popular QTFR is the 
Wigner-Ville distribution [13] for which the kernel is 
<K#,t) = 1: 



r +oa t t 

,/) = J z(t+-)z*(t--)e-^dr. (17) 



w z (t 



To reduce the artefacts arising from the cross terms 
between the different components of the signal, we 
use the analytic signal [2] z(t) defined by 



z(t) = x(t)+jH[x(t)} 



(18) 



where H[x(t)] = ^ * x(t) is the Hubert transform 
of x{t) designed by the Kaiser window method (FIR 
Hilbert transformer) [T5]. The Pseudo Wigner-Ville 
distribution is given by 



PW z (t,f) 



+ 00 



z{t+ T -)z*{t-\)e-^dTllV) 



where h(r) is a smoothing window (in our case a 
Kaiser function). 
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5 Experimental results 



5.2 Linear case 



5.1 Experimental apparatus 

Figure Q] shows the experimental setup. It consists 
of a 8 m long cylindrical pipe having a 5 cm inner 
diameter and a 0.5 cm thick wall connecting with an 
array of 60 Helmholtz resonators as side branches. 
The distance between two consecutive resonators is 
d n =0.1 m. The upstream section links this system 
to a loudspeaker designed for high acoustic power 
level and used to generate linear frequency modu- 
lated waves (chirps) or wavepackets (approximate 5- 
function). The duration of the shirp is 0.5 s and 
the frequency range is included in [100; 600] Hz. The 
wavepackets duration is 0.01 s and it frequency range 
extends from to 2500 Hz. At the end of the down- 
stream section, an approximately anechoic termina- 
tion made of plastic foam suppresses reflected waves. 
As noted above, the QTFRs are able to discriminate 
the reflected signal from the incident one and an ane- 
choic termination seems unnecessary. However, its 
use prevents a too heavy contamination of the use- 
ful picture in the time-frequency plane by unessen- 
tial signals and thus improves the accuracy of the 
QTFR. Lastly, two microphones mi and m 2 (B&K 
4136 with Nexus 2690 amplifier) measure the pres- 
sure in up and downstream sections. These micro- 
phones produce 0.2 % of distortion at 150 dB which 
ensures that the nonlinear effects are generated by 
the propagation under the lattice (and not by the 
microphones themself). The data acquisition is car- 
ried out by means of a 16 bits AD-converted with 
a sampling frequency of 10 kHz and an anti-aliasing 
filter (8 th order Chebyschev filter) with a bandwidth 
of 4 kHz . 

All the resonators are identical and are L = 16.5 
cm long cylindrical cavities with a diameter of 4 
cm. Their volumes may be independently tuned by 
moving pistons. Their neck is a 2 cm long tube 
having a diameter of 1 cm. The Helmholtz reso- 
nance frequency corresponding to the whole volume 
(V = 2.1 10~ 4 m 3 ) is 300 Hz. All the pipes used are 
sufficiently stiff that their structural modes are not 
in the frequency range of interest. 



Figure [2] gives the theoretical [/, cos (qd)] plot of the 
dispersion relation of Bloch waves for the ordered 
lattice described above (Fig. [2]-a) and the exper- 
imental results deduced from the transfer function 
of the system (Fig. [2|-b). We observe a good fit in 
the frequency range of interest. For this geometry 
the curve shows two kinds of gaps : those (labeled a 
and b) due to the resonators included in [300 : 450] 
Hz (Helmholtz resonance) and [1100 : 1200] Hz (cor- 
responding to the length of the cavity L) and one 
other that is a specific characteristic of the periodic 
lattice (dn = 0.1 m) at 1700 Hz (labeled c). The 
Pseudo Wigner-Ville Distribution images associated 
with the propagation of an acoustic pulse in the lat- 
tice are presented in Fig. [3) Figure[3]-a is the PWVD 
of the signal picked-up in the upstream section of the 
apparatus by microphone mi corresponding to the 
incident impulse beginning at 0.4 s and three narrow 
bands reflected signal. Fig. [3}b shows the PWVD of 
the transmitted (output) signal given by microphone 
m,2 (no reflection are visible because of the ancchoc 
termination). As well as the time frequency plot of 
the incident impulse, other details can be seen in Fig. 
131-a. For time t > 0.47 s a part of the input signal is 
reflected corresponding to the upper cutoff frequen- 
cies of the gaps. The two lower stopbands correspond 
to the energy of acoustic waves reflected by the res- 
onators. 

The third case has a different nature since this gap 
is due to the spatial periodicity of the lattice. The 
reflected signal corresponding to this stopband indi- 
cates that at least for a finite length lattice, the spec- 
trum structure (gaps and allowed bands) is not pro- 
duced only by wave interferences but also by some 
wave reflection. This result cannot be displayed by 
a method other than QTFRs since it requires time 
and frequency to be localized. Some authors have 
investigated analyticaly [3] or experimentaly [16] the 
dispersion characteristics of sound waves in a lattice 
of Helmholtz resonators but they have not pointed 
out the existence of a reflected wave for particular 
bands. 

In Fig. [3]-b, the waves belonging to the allowed bands 
are shown at the time 0.49 s corresponding to the 
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propagation of waves through the lattice. The gaps 
are easily identified and are in good agreement with 
the values given by the Bloch theory. The waves with 
frequencies closed to the gaps arrive at longer and 
longer delays : their group velocity goes to zero as 
the frequency gets closer to the cutoff frequencies. 
This is indicated on the PWVD plot by long tails on 
each side of the gaps. From the point coordinates of 
PWVD of the output signal, it is possible to estimate 
the group velocity as a frequency dependant function 

(Fig. rj. 



dB dB 




time (s) time (s) 



Figure 3: (a) Wigner-Ville transform of a pulse signal 
upstream of the lattice, (b) Wigner-Ville transform 
of pulse signal downstream of the lattice. 



5.3 Nonlinear case 

In the nonlinear case, the physical parameters of the 
output are dependent on the input parameters. It is 
then necessary to consider several cases with different 
input conditions. To investigate the effects of these 
nonlincarities, the input signal is a linear FM signal 
with a frequency range including the first stopband 
of the lattice (from 100 Hz to 600 Hz). Fig. [5] and 
Fig. [6] show the PWVD of this chirp for low and 
high acoustic levels corresponding to the linear and 
the nonlinear cases respectively. In the linear case, 
waves belonging to the first forbidden band cannot 
propagate, as shown by Fig. [SJ-b. For higher acoustic 
levels, nonlincarities generate harmonics 2u and 3a; 
which belong to allowed bands (Fig. ([7]) which repre- 
sents a zoom in the time scale range [0.2s : 0.4s] of the 
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Figure 4: Group velocity vs. frequency and transfer 
function of the lattice. 



Fig. ([6])). Then, for some particular frequencies in 
the first stopband the nonlinearities take part to the 
derealization of the energy and gives rise to energy 
transport : the energy of an input signal contained in 
a frequency range which is in a gap propagates in the 
lattice after being transferred in ranges correspond- 
ing to allowed bands in the spectrum. This mech- 
anism breaks the filter structure of the lattice and 
is partly responsible for changes in the transmissiv- 
ity of the medium. Indeed, the amplitude is restored 
with presence of nonlinearity : each site exhibits a 
nonlinearity and for localization lengths of the order 
of the lattice spacing, the associated wave becomes 
propagating. The measurement of the energy ampli- 
tude of each harmonics can permit the estimation of 
coefficients a n (and the next in the Taylor's develop- 
ment) corresponding to the harmonics in the model 
developed in section 3. The Fig. (|7|) shows also the 
different phenomena induce by the nonlinear effects : 
in the beginning of the stopband (for time including 
between 0.2 and 0.25 s) only harmonic 2 is present 
in the response of the lattice whereas in the rest of 
the stopband harmonics 3, 4 and 5 are detected in 
the time-frequency image. Because the nonlinear ef- 
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fects are not the same in all the stopband, we pre- 
sume that the different coefficients of the nonlineari- 
ties in the model developed in section 3 are depend- 
ing on the frequency of the wave : these effects result 
from the competition between the different orders of 
the nonlinearities [17 . A experimental study of one 
Helmholtz resonator confirms this assumption and al- 
lows the measurement of the different coefficients of 
the model [IS]. 
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Figure 5: (a) Wigner-Ville transform of a FM signal 
upstream of the lattice (linear case), (b) Wigner-Ville 
transform of a FM signal downstream of the lattice 
(linear case). 
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Figure 6: (a) Wigner-Ville transform of a FM signal 
upstream of the lattice (nonlinear case) . (b) Wigner- 
Ville transform of a FM signal downstream of the 
lattice (nonlinear case). 




Figure 7: Zoom of the figure ((6]b) in the time range 
[0.2 : 0.42] s. Numbers label the corresponding har- 
monics from the fundamental (1) to harmonic 5. 

6 Conclusion 

In this paper we have shown that the time-frequency 
methods are well adapted to the investigation of 
the propagation of dispersive waves. Especially they 
give information about the group velocity and allow 
to find the dispersion characteristics of complicated 
waveguide. These methods are also very efficient for 
displaying the nonlinear behavior of high-intensity 
sound waves : the generation of higher harmonics 
can be directly measured. Lastly, the interpretation 
of experimental results of wave propagation in non- 
linear lattices is made easier with the help of time- 
frequency methods. 
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